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ABSTRACT 


A digital data acquisition system was implemented for the study and collection 
of scattered waveforms from various simple objects. Reflected waveforms were 
extracted from interference with other signals by a simple subtraction technique. 
Backscattered and bistatic reflections were observed and recorded. Various effects 
on scattered waveforms, including target composition, pulse length, and frequency, 
were experimentally observed. Bistatic scattering results from a stryofoam sphere 
at a ka value of 10 were in agreement with theoretical predictions for a rigid sphere. 
Bistatically measured data at 70 kHz from a styrofoam cylinder of length 20 cm and 
diameter 2.5 cm followed a trend for bistatic angles up to 150 degrees but became 


amorphous as the bistatic angle approached 180 degrees. 
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I. INTRODUCTION 


Any object with a characteristic impedance different from the surrounding 
fluid will affect an acoustic wave impinging on it. The wave will be transmitted 
through the object and scattered from it. Flexural waves, compressional waves, and 
creeping waves are examples of some of the possible scattering and reflection 
mechanisms [Ref. 1:p. 301]. The “characteristic impedance’ equals the density of a 
material multiplied by the speed of sound in the material. It is one measure used to 
calculate reflection factors for pressure waves. Other properties that affect the 
scattering properties of an object include frequency, shape of the object, and pulse 
length. | 

Ever since the development of the submarine as an important force in naval 
warfare, much effort has been directed to improving their detection. The same 
holds true for the submarine detecting of surface ships or other submarines. The 
knowledge required for detection include knowing the source level of the target, 
propagation loss, noise level , and lastly target strength, the subject of this paper. 
These parameters are used in the well-known sonar equations [Ref. 2:p. 411]. 
Target strength, TS, as used in the sonar equations is defined as: 


J 
TS= 10 log— (1) 


O 


where Ir is the reflected intensity one meter from the acoustic center of the target 
and Ig is the incident intensity at the target without the presence of the target. This 
equation is general, giving no stipulation as to the direction of the reflected 
intensity. Typically, signals in active sonar are sent and recevied by the same 


platform, thus the most common geometry and ,therefore, the most studied, is the 


backscattered monostatic case where the reflected intensity 1s measured in the 
direction back toward the source. Measurement along any other direction 
constitutes a bistatic geometry. Relatively little effort has been given to the study of 
bistatic target strength until recently when the bistatic geometry was recognized as a 
potential method for improving underwater detection of submerged targets. 


Figures | and 2 illustrate the monostatic and bistatic geometries. 
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Figure 1. Backscattered Monostatic Geometry 
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Figure 2. Bistatic Geometry 


Some work has been done in the radar community on estimating the bistatic 
radar cross-section of objects. Crispin and Siegel [Ref. 3:p. 156] present a theorem 
for the approximation of radar bistatic cross-section: 

In the limit of vanishing wavelength, the bistatic cross-section for 


transmitter direction K and receiver fi is equal to the monostatic cross-section 
for the transmitter-receiver direction k + f with & # fi for bodies that are 


sufficiently smooth. 


(Figure 3 illustrates the geometry as stated in the theorem.) 
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Figure 3. Radar Bistatic Cross-Section Theorem Geometry 
Their bistatic theorem has its limitations when applied to complex shapes; The 
bistatic cross-section can be greater or less than the monostatic case due to 
structurally-related interferences. Another limitation is the size of the target with 


respect to the incident wavelength. Spheres were analyzed using various ka values 


where where a is the radius of the sphere and k is the wave number. It was found 
that the theorem held more closely for higher ka values, but in all cases the theorem 
sharply deviated from experiment as the bistatic angle approached 180 degrees. As 
an illustrative example, they presented data on a jet aircraft and a ballistic-missile 
shape [Ref. 3:p. 174-177]. Data were in fairly good agreement with their theorem 
as long as the bistatic angle did not get too large. 

In transition from the radar community to underwater acoustics, Urick [Ref. 
1:p. 302] implied that the bistatic radar cross-section theorem might apply in a 
underwater condition, in which case the bistatic target strength was equal to the 
monostatic target strength taken at the bisector of the bistatic angle. Urick 
emphasized there was a lack of data to prove or disprove such a theory. It would be 
very convenient if such a theorem existed since there is plentiful target strength data 
for monostatic geometries. Whether or not the high frequency limitations of such a 
theorem would preclude its tactical use is another question. 

The interest in bistatic target strength has been generated by the possibility of 
increasing detection of submerged and surface units using this technique. For 
example, one can easily visualize having a surface ship using its active sonar while a 
submerged submarine, using passive sonar, listens for reflections off enemy forces. 

There are several advantages and disadvantages to bistatic employment. Several 
ideas now discussed were taken from Harvey [Ref. 4] who gave a thorough 
discussion on using bistatic employment as a tactical tool in naval warfare. The 
major advantages include the ability to maintain a track of the enemy force using 
passive sonar and the increased amount of data available to detect and classify 
targets. All units besides the active unit can be passive, thereby hampering own- 


ship detection. Detecting bistatic reflections, as well as backscattered reflections, 


from a target increases the data available for detection and classification of the 
target. The greatest difficulty using a bistatic geometry is the high positional 
accuracy required between the pinger and the listener. This requires prior planning 
for use with surface ships and submarines since they do not normally have 
continuous open communication with each other. Another difficulty is the 
distinguishing of backscattered reflections from bistatic reflections. It may turn out 
that the advantages do not outweigh the disadvantages and the whole idea will prove 


infeasible. 


Il. RGET STRE H 


A. INTRODUCTION 
This chapter presents some definitions and background knowledge essential in 
preparation for experimentation. Some concepts are rather simple but are included 


for completeness. 


B. DEFINITION 

The definition of target strength has already been given, equation (1). 
However, there are several aspects of the definition that warrant clarification. 
First, all measurements are taken in the far field of the target. The far field is 
where the pressure reflected from an object exhibits a 1/r relationship [Ref. 2:p. 
178]. For example, pressure measurements taken one meter from the hull of a large 
ship would certainly be in the near field, so intensity is measured in the far field at 
ranges of hundreds of meters and then extrapolated back to one meter for use in the 
target strength calculations. 

The next question concerns the defintion of the acoustic center. Very rarely are 
we able to know the location of a target's acoustic center. The acoustic center is the 
point from which sound appears to radiate. It may be anywhere around the target, 
and may lie outside the shape itself. It may also change with the aspect of the target 
relative to the receiver. Exact location of the target's acoustic center is ignored in 
field measurements, and the geometrical center is used. When measuring pressure 
at hundreds of meters, a few meters error in determination accounted for by the 


acoustic center will not produce a large discrepancy. 


C. SPHERICAL WAVES 

For a distant source, to a good approximation, the wave will appear locally as a 
plane wave as it encounters the target. Consider an omnidirectional point source: 
Spherically diverging waves are emitted and the pressure amplitude and phase is 
dependent only on the distance r from the acoustic center of the source. An 


outwardly travelling spherical wave can be represented by 


A _j(t-k) 
@ 


P (r) (2) 


A = pressure amplitude at one meter 
r = distance from point source 

w = angular frequency of signal 

k = wave propagation number 


When the radius of curvature of the wave front is large compared to the biggest 
dimension of the target, the spherical wave can be treated as a plane wave. 

Plane waves are assurred when the dimensions of a target are small compared to 
the source-target distance. One of the goals of small-scale experiments is the 
application of the results to oceanic situations, therefore there is a concern about the 
approximation of plane waves ensonifying a target in a small laboratory tank. 

The geometrical constraints of the water tank limits target-source-receiver 
distances. With small distances between source and receiver and big targets, a plane 
wave may not be a good approximation. Another factor that influences the plane 
wave approximation is the directionality of the source. It 1s difficult to achieve 
omnidirectional behavior from a finite source and it is impossible if the frequency 


is high. Beam patterns are formed with energy emitted in preferred directions. If 


the target subtends an angle greater than the beamwidth, the pressure amplitude will 
not be the same for all points on the target. 

Plane waves also have properties that make target strength calculations simple 
to handle. The intensity can be equated to the root-mean-square pressure amplitude 
squared divided by the characteristic impedance [Ref. 2:p. 111]. Since hydrophones 
measure pressure via a Sensitivity factor, the output voltage they produce can be 


directly related to intensity. 


D. WAVEFORM EFFECTS 
There are several factors to consider in the measurement and collection of 
reflected waveforms from objects. Some are listed below and each will cause the 
reflected waveform to differ from the incident waveform. 
- target composition 
- target shape 


- pulse length 
- resonance effects 


For the reflected waveform to be identical of the incident waveform, a planar 
rigid boundary must be encountered. For more complicated targets, some of the 
incident energy will be transmitted into the body and eventually reradiated and 
different areas of the target may reflect in such a way to cause mutual interference. 
Many objects have resonances, at which the target will store some of the incident 
energy and reradiate it at a later time causing interference effects. Theoretical and 
experimental results show that for short pulse lengths, reflections are nearly 
independent of target composition and are almost identical to the incident pulse 
[Ref. 5,6,7] while for longer pulses the reflected waveform changes drastically 


depending on the composition and shape of the target. 


E. EXPERIMENTAL MEASUREMENT 

While conducting the literature search, it was interesting to note how the 
pressure amplitude of the reflected wave was measured. As Stated before, 
hydrophones give a voltage that is representative of pressure and this can be equated 
to intensity. The problem lies in the fact that the reflected waveform, in most cases, 
is not identical to that of the incident signal. The experimenter has to specify what is 
measured and what kind of pulse is being used. Response can be catagorized as 
transient or steady-state. Steady-state implies that the whole target has been 
ensonified by the incident pulse allowing for all portions of the target to contribute 
to the reflection and for any resonances to become fully excited. In this case, the 
return will be a combination of all possible reflection paths and modes. A transient 
response simply is not in the steady-state. To give an idea of measurement used, the 


following are quoted from two published articles: 


...five cycles of the sine wave located at a position in the pulse judged to 
be a true representation of the steady state. The average of the ten peak values 
so obtained gave the peak pressure amplitude [Ref. 7] 


After amplification, the amplitude of the back-scattered signal is 
measured with a peak detector [Ref. 8]. 

The main point is that for reflected waveforms, not similiar in shape to the 
incident waveform, some subjective determination is required. Estimation of the 
Steady-state must be made or some definite parameter has to be selected such as peak 
amplitude or total energy. In any case, an unique way to convert the amplitude 
measurements to an intensity level is not achievable. Add to this the fact that the 
reflected waveform will be different for different target orientations and 
frequencies, and the problem of translating measurements into an intensity level to 


calculate a target strength becomes complicated. 


There presently exists a Target Strength Standardization Program [Ref. 9] 
which was assigned the task of finding some standard measure of target strength 
independent of signal type. The general consensus from the program was to use the 
"Average Integrated Target Strength" which 1s the average time-integrated energy. 
Furthermore, an Echo Metrics Subcommittee [Ref. 9] was formed to find signal- 
processor dependent metrics for the measurement of target strength. Theoretically, 
the incident energy encountered by the object would be reflected or absorbed and 
eventually reradiated at some later time providing no energy was stored or lost in 
the object. Therefore, it would seem that the average of the integrated energy 
would be an independent metric for all types of signals providing the integration 
time is extended to include all the return signal. Some standardization of target 
strength is being investigated so that data from independent experiments can be 


compared and integrated. 
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A. BACKGROUND AND CONCEPT 

The experiments were conducted in a water tank located in Room O19 in 
Spanagel Hall at the Naval Postgraduate School. The tank was approximately 1.5 m 
wide, 8 m long, and 2 m deep and lined with anechoic rubber material. For the 
targets used, wall reflections were of the same magnitude as the target reflections. 
Usually, it was easier to use pulsed signals and time discriminate for the target 
reflection. However, the small dimensions of the tank did not provide the 
environment to do this effectively for all geometries. To overcome this difficulty, a 
subtraction technique was used employing a digital storage device. The basic 
concept was to subtract the signal received with the target immersed in the tank 
from the signal received with the target absent. The procedure was very simple and 


allowed for any target-source-receiver geometry. 


B. DEVELOPMENT 

A considerable amount of time went into the selection of the digital sampling 
and storage device, the source, and the receiver. The first question that had to be 
addressed was the range of frequencies to be used. The frequency range of 100-200 
kHz was initially selected because it gave wavelength-to-target dimension ratios that 
match oceanic scenarios. The critical specifications for the digital storage device 
are the sampling rate, dynamic range and jitter. Dynamic range was not a problem 
since most digital oscilloscopes have the same voltage range as a typical non-digital 
oscilloscope. The sampling rate must be fast enough to provide a good 


representation of the signal, and assuming a sinusoidal signal about ten samples per 


I] 


cycle would give a good trace. For example, fora 100 kHz signal the sampling rate 
would have to be 1 MHz. Another major concern was the amount of jitter in the 
sampling device. Jitter describes the fluctuation in the triggering times relative to 
the sampled signal. For a system using a subtraction process, it is perhaps the most 
critical specification because there has to exist a common time reference between 
collected waveforms to get accurate results. 

The selection of the digital oscilloscope was based on the availability of 
equipment. Two digital oscilloscopes were evaluated. The Nicolet 3091 was 
analyzed first. Its fastest sampling rate was 1 MHz and had a storage capability of 
4000 points. The Nicolet 3091 had a RS-232 serial interface and was easy to use, 
however it was not acceptable due to excessive jitter at the frequencies of interest. 
The jitter was easily seen by applying a pulsed signal and observing the traces on the 
scope. The digital oscilloscope finally used was the Tektronix 2430 (Tek 2430). It 
is a very capable osciiloscope and has many nice features making it suitable for data 
acquisition. The most important feature that it had was an averaging mode 
selectable from 2 to 256 waveforms. This helped eliminate random noise and 
minimized any fluctuations the electronic equipment possessed. The amount of 
jitter was dependent of the time scale selected. The fastest sampling rate of the Tek 
2430 was 100 MHz which could be used to sample a 10 MHz signal using ten points 
per cycle. The jitter specification given in the operating manual, for the time- 
delayed mode, was .08 x (time scale). The jitter was observed and was significant 
for the time scale desired, however the saving factor was the time averaging mode 
which ultimately averaged out the jitter and provided a sufficiently constant time 
reference. It was very easy to observe this on the oscilloscope because of a 


triggering marker on the signal. Jitter becomes less of a problem for any digital 
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oscilloscope when the number of samples per cycle of signal increases. However, 
the problem then arises from the limited amount of storage in the oscilloscope. A 
smaller number of cycles can be stored when using a faster sampling rate. The Tex 
2430 has a storage capabiltiy of 1024 points, much less that the Nicolet 3091; 
however it also has a time-delay feature which allowed for selective data storage 
(The Nicolet 3091 did not have any time-delay feature). Other features of the Tek 
2430 include cursor control for voltage and time measurements, simple output to a 
Hewlett Packard Think-Jet printer or pen plotter, and a GPIB interface. It is a very 
versatile piece of equipment and was finally selected for use in the data acquisition 
system. 

The selection of the source and receiving hydrophone was also limited to in- 
house equipment. For the frequencies desired, the USRD type F41 transducer was 
the appropriate choice with a nominal frequency range of 15-150 kHz [Ref. 10:p. 5]. 
It is a piston type piezeoelectric transducer with a face diameter of 10 cm. The 
choice of the receiver was a Celesco LC-10 hydrophone. It has a wide nominal 
frequency band, is small and light making it easy to move around, and has good 


sensitivity. 


C. SYSTEM SET-UP 
1. Electronic Equipment 
The F41 source had to be driven near 200 V rms for target reflections to be 
seen. An Instruments, Inc. Model L2 power amplifier was used to drive the source 
and its selection was again based on what was available in the laboratory. It allowed 
for impedance matching and has a voltage-divided 10 mV/V port to monitor output 
voltage. The output of the receiver was passed through an Ithaco Model 1201 Low 


Noise preamplifier which has bandwidth setting capability. This was very 
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beneficial because there was a lot of 60 Hz noise which had to be eliminated by 
bandwidth filtering. The signal out of the preamplifier was sent to a Kikusui 
COS5060A oscilloscope and the Tektronix 2430 digital oscilloscope. Having a non- 
digital oscilloscope made observing and measurement of the waveforms easier to 
perform on the Tek 2430. A Hewlett Packard 3314A function generator was used 
to provide the source signal. It also provided the triggering signal which was sent to 
the external channels of the oscilloscopes. In the pulse mode, the HP 3314A sends 
out a3 V transistor-transistor logic signal via its trigger port. The pulse repetition 
rate of 10 ms was sufficient to allow reflections to diminish before the next pulse. A 
Hewlett Packard 9000 Series 300 desk top computer was used as the data collection 
and processing unit. It was connected to the Tek 2430 via the HPIB bus. The 
transfer of waveforms stored on the Tek 2430 was controlled by software 
programs. Data was ASCII encoded from the Tek 2430, transferred, decoded into 
actual voltages, and then stored onto disks in binary data form. A schematic 
diagram of the electronic system is shown in Figure 4. 
2. Physical Set-up 

To measure all possible bistatic angles, the scattering pattern would have to 
obtained in three dimensions. This was not done in this work. Bistatic data were 
limited to the horizontal plane determined by the source, receiver, and target, with 
the axis of the target in the same plane. 

The source, the F41 transducer, was mounted on a 3/4-inch aluminum rod 
attached to a rectangular support plate. The rod was fixed to the plate by means of a 
friction fitting secured by a hose clamp. Target, source, and receiver were 
submerged one meter below the surface of the water and the rods which supported 


them were marked for this depth. The source with its support rested on another 
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Figure 4. Data Acquisition System 
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support built from aluminum L-beams that spanned the width of the tank. Lead 
bricks were used to hold the source once it was set in place and a level was used to 
ensure the source face was vertical. The receiver was mounted onto a sliding track 
specially built for the tanks. Adjusting screws were added for ease of positioning 
the receiver vertically. The initial intent was to design a rotating boom to hold the 
receiver and support the target from its center. A device was constructed, but too 
much receiver movement occurred between the loading and unloading of targets. 
Though the boom was never used, the device was still used for the suspension of 
targets and had a mechanical angle indicator on its top for angle measurement. 
Spherical targets were suspended by a thin cable and the cylindrical target by a 3/8- 
inch rod. 

The direction of the source was Set by aligning the source, receiver, and 
target along a straight line with the receiver between the source and target. Source, 
receiver, and target depth were one meter. Vertical alignment was ensured by 
keeping the source and receiver vertical by use of a level. The frequency was set to 
200 kHz to establish a narrow major beam and the source was rotated until 
maximum signal was seen by the receiver. This ensured the target was being totally 
ensonified by the major lobe of the source. 

Once the source direction was set, the receiver was moved to the desired 
position. The angle between the source, target, and receiver was determined by 
setting the receiver's track perpendicular to the source-target axis and then using 
trigonometry relationships. After the desired geometry was Set, care was taken to 
ensure the target, source, and receiver were stationary. It was critical that all 


elements be stationary when waveforms were collected. 
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3. Data Collection 


The method of data collection and storage followed the sequence: 


set up desired source-receiver-target geometry 


warmup all electronic equipment >1 hr 


ensure connections as illustrated in Figure 5 
load "Tex" onto computer 
set up the Textronix 2430 oscilloscope 


- B trigger mode to "runs after" 
- A trigger 

- mode: normal 

- coupling: AC 

- source: external 1 

- slope: positive 

- level: 800 mV 

- trigger position: 1/8 
data encoding: ASCII 
GPIB address: 4 

vertical mode: channel 2 


select tume delay to obtain desired time window 
set acquisition mode to 256 averaging 


adjust channel 2 voltage scale/Ithaco gain appropriately 
for the signal 
for each new signal, press the acquire button 


wait sufficient time for averaging to take place 


press the save button 
run the "Tex" program 


The "Tex" program controlled the transfer and storage of data from the 
digital oscilloscope. Inputs from the user were file names and file descriptions. 
The transfer process, including writing onto the disk, took approximately 10 s per 


waveform. A listing of the program is in Appendix A. 
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The computer program "Plot" was used to evaluate the stored data. It 
took two files and subtracted them, and then plotted the results on the CRT 
terminal. Hard copy output could be obtained by porting to a pen plotter or Think- 
Jet printer. Input from the user required file names and where printout was to 
occur. The plotting routine was written so that the plot's parameters were easily 
changeable. A listing is enclosed in Appendix B. 

The method of data collection was a slow and somewhat tedious process 
since all movement of target and receiver ware manually controlled. Data 
readings were also recorded manually. 

4. Validation 

The question of whether the data acquisition system could correctly 
extract a target reflection immersed within other signals had to be answered. The 
very simple procedure of subtracting background noise works only if there is a 
common time reference for "target in’ and "target out" waveforms. For ease of 
discussion, “background noise” is used for any signal, random or coherent, that is 
not the desired signal and corresponds to the ‘target out" condition. 

The procedure to validate the system began with obtaining a geometry 
which allowed a direct pulse and reflected pulse to be completely separated by 
time discrimation with no interfering wall reflections. Figure 5 is such an example. 
The larger pulse was the direct path pulse and the second pulse is the pulse 
reflected from a styrofoam sphere. The frequency was 66.5 kHz and the pulse 
length was 20 cycles. The reflected pulse is shown enlarged in Figure 6. The 
peak-to-peak voltage of the reflected pulse was measured to be 68 + 2 mV using 


the cursor controls of the Tek 2430. 
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The pulse length was then increased to 120 cycles so the direct pulse would 
be interfering with the reflected pulse. No other changes were made. The 
interference of the two signals is shown in Figure 7 and was stored in file A. Next, 
the sphere was removed and the background noise sampled and recorded in file B. 
The background noise in this case was mostly the direct pulse. The last step was to 
subtract file B from file A to get the reflected signal from the sphere. The result is 
shown in figure 8. The waveform begins at the same time in comparison with 
Figure 6 and the only difference is that the pulse length was longer as expected and 
the curve slightly more choppy in form. Measured voltage was 73 +4 mV. The 
larger uncertainty arose due to the combination of small jitter and amplitude 
fluctuations of the HP 3314A. The uncertainty was estimated beforehand by taking 
identical signals directly from the function generator to the Tek 2430, recording 
them in separate files, and subtracting them and observing the results. It was found 
that peak-to-peak voltage uncertainty was approximately four digitizing levels 
regardless of the voltage scale used (a digitizing level equals the voltage per division 
divided by 50 for the Tek 2430). Given that this error would always be present, 
there was good agreement between Figure 6 and Figure 8, and the method of 
obtaining data was considered validated. 

As with any signal processing system, a high signal-to-noise ratio is 
desired. For any digital subtraction routine, the desired signal has to be large 
enough relative to the noise so it will not disappear in the digitizing level of the 
sampling device. For example, if noise plus the signal were 10 volts using 512 
digitizing levels or bins, each bin represents 19.5 mV. If the signal were small, 1.e., 
1 mV, the signal level would be much smaller that the digitizing level and would be 


lost in the subtraction routine. 
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Figure 5. Direct and Reflected Pulse 
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Figure 6. Reflected Pulse 
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Keeping the source and receiver positions fixed relative to each other was 
critical. Small movement between the two greatly affected the phase of the received 
signal and had a large effect on the results. This effect was verified in preliminary 
work. The initial intent was to place the receiver and support the target on the same 
platform. However, the removal of the target for the background noise 
measurement caused flexing in the platform and produced source-to-receiver 
movement giving incorrect results in the subtraction routine. Figure 9 illustrates 
the effect of source-receiver movement on “target-in" and "target-out" waveforms. 
The solid trace represents the target-in waveform and the dashed line represented 
the target-out waveform. The reflection began at 35 usec. The dashed and dotted 
line is the the result from the subtraction of the two waveforms and it was clear that 
in this case, it did not give the correct results. The region between O and 35 usec 
should have cancelled to zero. The non-zero result indicates a phase change caused 
by relative source-receiver movement. The phase change subsequently affected the 
amplitude of the reflected signal from 35 to 100 usec. It was easy to see that any 
phase change brought about by small source-receiver movement would cause 
erroneous results. To eliminate this problem, a separate platform was subsequently 
built and used for the receiver. 

The output of the power amplifier was also checked in preliminary work 
for fluctuations. An arbitrary source-receiver geometry was selected without any 
target in the tank. Several waveforms of the direct pulse were stored at various 
times after initial start up. Amplitude fluctuations existed shorty after start-up. 
Ideally, successive waveforms subtracted from each other would cancel out to give a 
zero response. If amplitude fluctuations existed a false signal would occur and this 


in fact did occur after five minutes of warm-up as Figure 10 illustrates. Another 
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test taken at 15 minutes after start up produced better but still unsatisfactory results. 


After one hour of warm-up the amplitude fluctuations subsided and a zero respose 


was obtained. 
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Figure 10. Waveform Amplitude Fluctuation Effect After 10 Minutes 
of Warm-up 


IV. FAR FIELD DETERMINATION 


The measurement of the far field of several objects was conducted. An equation 
that gives the approximate beginning of the far field is given by [Ref. 2:p. 188]: 
2 
e 


— 3 
Tf = (3) 


L = the longest dimension of object 
dX. = wavelength 


It is always best to determine the far field by measuring pressure amplitude as a 
function of distance from the object. The equation for a diverging spherical wave, 
as described by equation 2, can be manipulated by taking the log of both sides to 
transform it into a linear equation with a slope of minus one. The far field can be 
determined experimentally when the observed slope of the curve starts to deviate 
from minus one. 

The far field of the F41 was determined to begin at approximately 100 cm and 
70 cm for 176 kHz and 70 kHz, respectively. The far field for the sound reflected 
from the styrofoam cylinder positioned broadside to the source was determined to 
begin at 40 cm for 70 kHz. The three-inch diameter styrofoam sphere's far field 
began at less than 20 cm at 70 kHz. Figures 11 through 14 are graphs of the data 
used to determine their far fields. Lines of slope equal to minus one are drawn on 


the graphs. 
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Figure 11. Receiver Voltage vs Distance, F41, 176 kHz 
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Figure 12. Receiver Voltage vs Distance, F41, 70 kilz 
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Figure 13. Receiver Voltage vs Distance, Cylinder, 70 kHz 
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V. SPHERES 


A. INTRODUCTION 

To achieve an understanding of the complexities involved in analyzing 
scattering effects of various objects, it was instructive to begin with a very simple 
and symmetrical object, a sphere. Several solid metal spheres were used. Spheres 
were aluminum and brass, but, the exact composition was not known or investigated 
since they were to serve only as preliminary tools in studying the general properties 
of waveform scattering. 

Since metal spheres are not perfectly ngid, they are capable of absorbing 
energy and subsequently reradiating it. Also, they possessed resonances. Many 
experiments in the past 30 years have been conducted with spheres and long 
cylinders, and the Resonance Scattering Theorem (RST) has evolved from this work 
[Ref. 5,6,7,8,11,12,13,14]. This theorem predicts scattering patterns which entail 
both the rigid-like properties and the resonance properties of an object. It was not 
the intent to examine this theorem in detail but only to observe some of its 


predictions experimentally. 


B. MONOSTATIC FORM FUNCTION 

The goal of replicating some of the data in the literature was sought to establish 
confidence in the experimentation process. Neubauer, Vogt, and Dragonette [Ref. 
7] conducted an analysis on elastic spheres using steady-state signals and observed 
the backscattered signal form function from an aluminum sphere as a function of ka. 


The reflection form function is defined as: 


p20 
a 


Ps 
Po 


(4) 








Ps = scattered pressure amplitude measured at r 
Po = incident pressure amplitude at sphere center 
a = sphere radius 


Their experimental data were plotted and compared against theoretical results 
and good agreement occurred. 

Part of Neubauer's experiment was repeated. A four-inch diameter aluminum 
sphere was ensonified by 50 cycle pulses from 100 to 130 kHz, and the backscattered 
pressure amplitude was measured. Time discrimination of the reflection was used. 
The amplitude of the reflected pulse was taken at what was judged to be the steady- 
State part of the signal. The receiver was 1.60 m from the center of the sphere, 
placing the receiver well within the far field of the sphere. 

The wavelength used in the calcuation of ka was dependent on the speed of 
sound in the water, which in turn was affected by temperature, pressure, and 
Salinity. For the fresh water tanks, salinity was ignored. The speed of sound was 
based on an empirical formula that applied to distilled water [Ref. 1:p. 107]. 
Although the water was normal tap water, the distilled water formula was a close 
enough approximation [Ref. 15]. Water temperature of the tank was a consistent 17 
degrees Celsius throughout the entire experimentation phase. At this temperature 
and normal atmospheric pressure, the speed of sound was calculated to be 1473 m/s. 
Experimentally, two LC-10's were placed on a track 180 cm apart and the 
successive arrival times of the pulse were measured via the Tek 2430 with the 
cursor controls. The speed of sound thus was measured to be 1480 + 20 m/s, placing 
1473 m/s well within the experimental uncertainty, and this value was used as the 


speed of sound in water. 


gy). 


The shadowing effect of the LC-10 in front of the F41 source was taken into 
account by placing a dummy LC-10 at that position and recording incident pressure 
at the sphere center location. Additionally, since the LC-10 was not perfectly 
omnidirectional, the support rod was marked to maintain a constant receiver aspect 
with the target. Placing the LC-10 in front of the F41 also caused additional 
reflections to occur. The direct pulse was reflected by the LC-10 back to the 
source's face and again reflected there back to the LC-10. More "noise" was present 
when measuring the monostatic case but appropriate location of the receiver 
eliminated any interference from those reflections. 

The data matched Neubauer's results well, and the form function as a function 
of ka is plotted in Figure 15. Data from Neubauer's experiment is also plotted in 
Figure 15. The only difference between results was a small shift in ka of 
approximately 0.2-0.3, not surprising since the compostion of the sphere was 
unknown. The minima corresponds to resonant frequencies of the sphere, where 
the specularly reflected wave is interfering with the resonantly radiated wave from 


within the body. 


C. WAVEFORM EFFECTS 

The effect of a sphere's resonance was examined in more detail. A four-inch 
diamter aluminum sphere was ensonified by a 30 cycle pulse at, above, and below a 
resonant frequency of 82.7 kHz. At resonance, the ringing of the sphere is clearly 
seen in Figure 16 when the direct pulse ends. Pulse length was approximately 367 
Us long. The sphere's stored energy is simply being reradiated. Figures 17 and 18 
show the waveform response 2 kHz above and below resonance. They are 


decisively different and clearly illustrate the resonance effect. 
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The effect of sphere composition on reflected waveforms was observed. A 
four-inch diameter brass sphere was ensonified by the same frequencies as the 
aluminum sphere. Brass has a characteristic impedance approximately 2.5 times 
that of aluminum. Figures 19, 20, and 21 show the reflected waveforms having 
frequency dependence that are markedly different from the aluminum sphere. 

The effect of pulse length was observed using a 4-inch aluminum sphere at and 
off a resonant frequency. Figure 22, at resonance, and Figure 23, off resonance, 
show reflected wavforms using a short pulse of 2 cycles. They appear very similar 
in shape. Figures 24 and 25 show the results for a pulse length of 5 cycles. Again, 
the results of the two are similar. When pulse length reaches 15 cycles, the effects 
of the resonance begins to be clearly observed, as Figure 26 and 27 illustrate. 

For comparison, the same procedure above was performed on a four-inch 
diameter brass sphere. Results are illustrated in Figures 28 through 33. The brass 
sphere did not show any resonance effect as determined in the preceding 
composition experiment for these frequencies, so its results were very similar to the 
aluminum sphere's results off resonance. Waveforms were practically identical for 
the shorter pulse lengths regardless of frequency and composition implying the 
reflections are purely specular. The results clearly demonstrated the dependence of 


pulse length on reflected waveforms. 


D. SELECTION OF STYROFOAM 

Since defining target strength becomes difficult when using elastic spheres, 
which was not the intent of this work, the need existed to do one of two things: (1) 
work at a frequency that gives good replication of the input pulse, or (2) find a 


target that gave a very high impedance mismatch. 
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Figure 16. Aluminum Sphere Reflection at Resonance, 82.7 kHz 
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Figure 17. Aluminum Sphere Reflection, 84.7 kHz 
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Figure 18. Aluminum Sphere Reflection, 80.7 kHz 
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Figure 19. Brass Sphere Reflection, 80.7 kHz 





Figure 20. Brass Sphere Reflection, 82.7 kHz 





Figure 21. Brass Sphere Reflection, 84.7 kHz 
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Figure 22. Aluminum Sphere Reflection, 82.7 kHz 
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Figure 23. Aluminum Sphere Reflection, 84.7 kHz 
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Figure 24. Aluminum Sphere Reflection, 82.7 kHz 
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Figure 25. Aluminum Sphere Reflection, 84.7 kHz 
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Figure 26. Aluminum Sphere Reflection, 82.7 kHz 
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Figure 27. Aluminum Sphere Reflection, 84.7 kHz 
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Figure 28. Brass Sphere Reflection, 82.7 kHz 
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Figure 29. Brass Sphere Reflection, 84.7 kHz 
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Figure 30. Brass Sphere Reflection, 82.7 kHz 
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Figure 31. Brass Sphere Reflection, 84.7 kHz 
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Figure 32. Brass Sphere Reflection, 82.7 kHz 


5 is 


N 


be my (i cca ah ey (mara) ie em in od Eo ae pea Poa go 


os 


Come 
—> 


t 
it 


ee 
———— 
——— 


-— 
0 ee ee 


ew ee ee 


we 
TE eee 
—_> 
—— er ee 
oF _ . 
——= 5 amt, 
— —* —_— 
= _. 
— 
— = 
ae ee 
aan. 
——— 
i 
i 
——EEE 
Woes c=, = 
ay OO cere 
ene 
Se 
i ge eg ce ng ge = 
a eae 3 rs ees 
_— 
— —— 
—— -— — sae Se 
eee 
ee 
ON, oe A ee 
oo 
= ———_ foe 


+r 


OF eS SE 


- — 
a 
2 op Se agp nr We reg 
es ee oe > 
—_ ci 


a —— 


NE TS nae. 


246 


Cmicrasecands ! 
Pie iene | iene 


u! 
cy) 


ea). 5 


Figure 33. Brass Sphere Reflection, 84.7 kHz 
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The most important factor in achieving a good reflection is the impedance 
mismatch between the fluid and object. Expemmental work done with metal spheres 
by Hamptom [Ref. 5] shows the effect of compostion using several metals. A 
pressure wave in water will see little resistance to a hard rubber sphere because its 
characteristic impedance is near that of the water. An object that gives a high 
impedance mismatch is an air bubble. The impedance ratio of water to air is 
roughly 3000. For the case of an air bubble, the boundary condition between the 
water and air is termed "pressure release’, since the pressure transmitted into the air 
is nearly zero [Ref. 2:p. 127]. The difference between a pressure-release sphere and 
a rigid sphere is best seen graphically in Figure 34 [adapted from Ref. 15: p. 359]. 
In the ka region used in the experiments, both the rigid and pressure-release form 


functions will be nearly identical and equal to one. 


form function 





2 A 6 8 10 12 


Figure 34. Backscattered Form Function of a Rigid and Pressure 
Release Sphere 


A styrofoam ball was used to approximate an air bubble not near resonance. 


The reflected waveforms were generally identical to the input pulse at all 


frequencies. This made styrofoam an ideal target for this work. Figures 35 
through 40 show reflections from a three-inch diameter styrofoam sphere using 


various pulse lengths at 100 and 106 kHz. 


E. BISTATIC SCATTERING 

A 360 degree pattern of the pressure amplitude scattered from a styrofoam 
sphere was obtained to compare the results with theoretical predictions for a 
perfectly reflecting sphere. Since a sphere is perfectly symmetrical, it was not 
necessary to take the entire 360 degrees of data but only 180 degrees of it. 

Theoretical equations have been derived and published by several authors that 
predict scattering of acoustic plane waves from a sphere [Ref. 16,17]. They tend to 
be quite complicated functions involving sums of Hankel functions. Figure 41 was 
adapted from Medwin and Clay [Ref. 16:pp. 192-193] and shows computed 
scattering patterns in relative amplitude for a perfectly reflecting sphere in water 
for several ka values. The incident plane wave is arriving from the nght. The 
patterns show that for increasing ka the major forward lobe gets larger in amplitude 
and narrower in width. More structure, i.e., secondary lobes, develop as ka 
increases. 

Bistatic data were obtained for a ka =10 using a three-inch diameter styrofoam 
sphere and are shown in Figure 42. The incident wave arrives from the bottom. 
Frequency was 61.5 kHz and an eight cycle pulse was used. The graphed data 
revealed some structure in the forward scattering region and a large forward 
scattering lobe. Comparison with the theoretical prediction for a ngid sphere 
(dashed line in Figure 42) show that the salient features of the scattered pattern 


were duplicated. 
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Figure 35. Styrofoam Sphere Reflection, 100 kHz 
N=5 


en —_ i. ——a 
+ 
f, ‘ f f : 
aininini : 
ae i‘ ene 
a i). k5k} ~ ic, | | a } | | | | | J re ie si as av pee 8 et OG eer ee 
= ne Mh a ky ‘ v , 
Sa lib ata aoa id ae -| 
1 i : (} f {! 4 
i i if l 
en 
-, Ge 





Te 3& i6ea 
(microsacands } 


Relative Time 


Figure 36. Styrofoam Sphere Reflection, 106 kHz 
N=5 
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Figure 37. Styrofoam Sphere Reflection, 100 kHz 
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Figure 38. Styrofoam Sphere Reflection, 106 kHz 
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Figure 39. Stvrofoam Sphere Reflection, 100 kHz 
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Figure 41. Theoretical Scattering Patterns for 
Rigid Spheres in Water 
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Figure 42. Bistatic Form Function of a Styrofoam Sphere 
ka = 10 
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VI. BISTATIC SCATTERI OF A STYROFOAM INDER 


A styrofoam cylinder of diameter 2.5 cm and length 20 cm with flat ends was 
used for bistatic scattering measurements. The particular diameter was chosen 
because the styrofoam was commercially available in that size. The length was 
selected using two criteria: (1) the length was small enough so that the receiver 
could be placed in the far field of the sound reflected from the cylinder and (2) it 
was large enough to clearly present a cylindrical shape. 

Since the F41 was not omnidirectional, there was concern of approximating a 
plane wave at the target. The lower the frequency, the better the pressure 
distribution across the target because of the larger major lobe, but the wavelength- 
to-target length ratio would be less and the signal-to-noise ratio also became poorer 
at lower frequencies. A frequency of 70 kHz was chosen at a source-to-target 
distance of 125 cm. Figure 43 shows the incident pressure distribution as a function 
of distance off the axis at a distance of 125 cm. The lower curve is for the 
monostatic case where the receiver shadowed some of the incident energy with a 
receiver-to-target center distance of 50 cm. A two percent pressure amplitude 
difference existed for the shadowed case. An eight percent decrease in pressure 
amplitude existed when not shadowed. 

The cylinder was mounted to a 1/8-inch diameter metal tube which had a 2 X 
1/2 X 1/16-inch flat plate epoxied to it to provide a means of rotating the cylinder. 
The flat plate was forced into the interior of the cylinder and provided a snug fit. 
The effect of the plate on the reflections was assumed to be neglible since it was 
completely surrounded by the styrofoam. The tube was as small as possible to 


minimize reflections. A string with weights was attached to keep the styrofoam 
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Figure 43. Pressure Distribution at 70 Khz 
cylinder submerged. The thin tube was soldered into a 3/8-nch diameter tube 
which supported the entire arrangement. The cylinder was level before entering the 
water and appeared to remain that way when in the water. Figure 44 illustrates the 
Support setup . 

Data were recorded for nine bistatic angles. Time discrimination of the 
reflected pulse was used whenever possible. The source-to-receiver distance was 
125 + 2 cm in all cases with the F41 being driven at 177 V rms. The Ithaco 
preamplifier setting was 20X with a 3-300 kHz bandwidth. Receiver-to-cylinder 
center distance varied, but the receiver was always in the far field of the cylinder. 
All data were extrapolated to one meter from the target center for comparison. A 
20 cycle pulse was used. 

The cylinder was initially aligned broadside to the source by eye to within an 
uncertainty of about + 4 degrees. Assuming symmetry, data were taken for 180 


degrees of rotation and several points were taken beyond that limit. Figure 45 
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Figure 44. Support Platform Setup 
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illustrates the angles measured in the experiments. The varied angle, A, is the angle 
made between the cylinder's normal and the line connecting the center of the target 


to the source axis. Angle B is the bistatic angle. 
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Figure 45. Measured Angles in Experiment 

The results of the experiment are shown in Figures 46 to 55. In most cases, 
strong peaks were observed when angle A was one-half the bistatic angle and 
secondary peaks with smaller amplitude occurred 90 degrees from that. This 
geometry seems to indicate specular reflection off the long cylinder's surface and its 
flat ends respectively. As the bistatic angle increases the secondary peak broadens 
and eventually disappears. Well into the forward scattering region, greater than 
150 degrees of bistatic angle, the data became less predictable The distinction of the 
primary and secondary peaks became less clear with additional peaks showing up. 
At a bistatic angle of 180 degrees no distinctive peaks were present and amplitude 


became very large in comparison to other bistatic angles. 
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Vil. CONCLUSION AND RECOMMENDATIONS 


A data acquisition system using the Tektronix 2430 digital oscilloscope was 
implemented to extract scattered waveforms in a noise-filled water environment. 
The system was verified to work and reflected waveforms were collected. 
Preliminary work revealed that the relative movement of the source and receiver 
between collection of waveforms was critical in obtaining accurate data using a 
digital sampling system with a subtraction routine. Data collection was slow and 
tedious. 

Several aspects influencing scattered waveforms, including target composition, 
frequency, and pulse length were observed and recorded. The effects of these on 
the waveforms of the reflected pulse produce difficutlies in measuring the intensity 
of pulses with waveforms not identical of the incident pulse. 

Styrofoam targets were used to circumvent the problem of waveform 
dependence on frequency, pulse length, and composition. Bistatic measurements 
were taken using a sphere and a cylinder as targets. The data from the sphere 
compared well with theoretical predictions for a perfectly reflecting sphere. 
Readily available data were not found for comparison with the results obtained for 
the cylinder. The goal of matching wavelength-to-target length ratios observed in 
oceanic conditions was not done in this work. The USRD type E27 transducer 
would be an excellent source for this application since it has a frequency range 100 
to SOO kHz. Its small size makes its beamwidths very large which is beneficial when 
receiver-target distance is limited by tank dimensions. The E27 is also very 
lightweight and easily moved. The effect of the jitter would have to be evaluated at 


the higher frequencies. 
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The bistatic radar cross-section theorem could not be compared to the data 
collected in this report. The wavelengths used in the experiments are not close to 
the "vanishing limit" as stated in the theorem. One common trend between the 
bistatic data and the theorem was that the data became less predictable as the bistatic 
angle approached 180 degrees. 

For future work, a physical structure must be made to facilitate positioning of 
the source, target and receiver in three dimensions so that data collection can be 
completely automated and accuracy improved. The rotation of supports should 
have a digital device to allow control by the computer. The design and construction 
of such a system would not be a trivial project and would require a great amount of 
ingenuity and experience with electronic devices. The payoff of having such a 
system would be great to students studying underwater acoustics and in particular to 
anyone studying scattering effects. Once a positioning system is built and 
completely automated, data collection will become simpler and more time can be 


devoted to analyzing it. 
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APPENDIX A 
COMPUTER LISTING OF "TEX" 


| SHARE TR HE HH et 
{ 


!'PROGRAM TEX 
| 


DATA AQUISITION FOR BI-STATIC TARGET STRENGTHS 
USING THE TEXTRONIX 243@ DIGITAL OSCILLOSCOPE 


! 

| 

! 

| 

| THIS PROGRAM TRANSFERS WAVEFORM DATA FROM THE TEX 2430 
! TO THE COMPUTER WHERE IT IS CONVERTED TO ACTUAL VOLTAGE 
| VALUES AND STORED ON A FLEXIBLE DISK. INPUT FROM THE 

! USER REQUIRES FILE NAME, DESCRIPTION, AND CONTINUATION 
! RESPONSES. PREVIOUS SETUP OF THE TEX 2430 IS REQUIRED. 
! A THOROUGH KNOWLEDGE OF THE GPIB CODES IS NECESSARY. SEE 
| THE OPERATOR'S MANUAL AND THE INTERFACING MANUAL FOR 

| THE TEX 2430. 

| 


| HHH HRI 
| INITIALIZATION 

| SHAR SHRI HEIR HHH RIES 
PtNreR IS CRI 

OPTION BASE @ 
DIM Vdata$(15] ,WS[ 40]! NAME OF DATA FILE AND DESCRIPTION 

DIM AS[ 200} ,BS[ 5000} ,DSI 5] !PRE-WAVEFORM AND WAVEFORM STRINGS 
DIM V( 1023) ,Tim( 1023) ,Volt( 1823)!DATA ARRAYS 

DIM X${ 10} ,¥Ys$[ 10} ,Vs$[ 10] , Ts$[ 12} 

DIM An$[ 1] 

| 
Bicsctccoccsc:s:ss:s:et:ssissers:tsiectrsraisssiastes:ereatsrers:eeseiatsiesasssrs:ssdte:eseisiaretscsterarsesessat 
TEXTRONIX SETUP 

| SPHERE HHH i ett 
; 


OUTPUT 704;"PATH OFF"!SHORTENS THE DATA STRING SENT BY THE 2430 
TST RE eS Hat et 
| INPUT FROM USER 

| HRUHHHHHHRHHRAHHHHHHHEHHHEHHHHHEH IHRE 
; 

INPUT “NAME OF THE DATA FILE?” ,Vdatas 

INPUT “DOES FILE ALREADY EXIST?¢(Y=YES)" ,Ang 

IF An$S="Y" THEN PURGE Vdata$ 

INPUT “INPUT FILE DESCRIPTION?(4@ CHARACTERS MAX)" WS 

BEEP 

PRINT "ALLOWING DATA DISK INSERTION. PRESS CONTINUE TO RESUME.” 
PAUSE 
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500 
510 
520 
53Q 
540 
55@ 
560 
57 
580 
59@ 
5OO 
519 
520 
530 
5b4@ 
65@ 
6560 
67@ 
58Q 
69@ 
70 
710 
72@ 
730 
7AQ 
750 
76@ 
77 
780 
79@ 
BOO 
81@ 
87Q 
8 3@ 
B4@ 
85@ 
86Q 
87@ 
88Q 
8390 
900 
41 
92Q 
3 3Q 
940 
35 
36@ 
97@ 
38Q 
I3@ 
1 OOD 


CREATE BDAT VdataS,1,17000 

ASSIGN @Dfile TO Vdata® 

PRINT 

IMHRHHHHHHAHHHHAHHHRRHRAHHHH HHH 
| OBTAIN WAVEFORM 

| SHH A HAR Rae a et inti 
! 

PRINT “OBTAIN DESIRED WAVEFORM ONTO TEXTRONI X" 

PRINT "PRESS CONTINUE TO RESUME" 

PAUSE 

1ST REE ot et aa 
| RECEIVING PRE-WAVEFORM DATA AND THE WAVEFORM 
HHS Ra at 
! 

OUTPUT 7@4;"WFM? XIN"! TIME INCREMENT BETWEEN POINTS 

ENTER 704; X% 

OUTPUT 704;"“WFM? YMU"!IVOLTAGE DIGITIZING INCREMENT 

ENTER 704;Ys$ 

OUTPUT 704;"“CH1? VOL" !VOLTS/DIV 

ENTER 704;Vs$ 

OUTPUT 704;“HOR? BSE"!B TME.SCAEE 

ENTER 704;TsS 

| 

OUTPUT 704;"CURV?" 

ENTER 704;BS 

| FHA eee eat 
SEPARATION OF ASCII STRING INTO DATA POINTS 

| HHHHRRHHHHHHHHHHHHHHHHHHRENHRHEHHHHHHHHAH HHH tas 
! 

N=LEN( BS) 

D$="" 

Inc=@ 

FOR I=1 TON 

LF BSlLIFiiS teen 

V( Inc)=VAL( DS) 

Inc=Inct1 

D$=""" 

GOTO 95@ . 

END IF 

DS=DS&BSII 1) 

NEXT I 

| SSH SR RRR et at 
| CONVERSION TO REAL-TIME DATA 


| SESeSeRee ett St TR at 
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101@ 
1029 
1@3@ 
104@ 
1@S0 
1@62 
197@ 
198@ 
1@9@ 
1100 
111 
112@ 
113 
114@ 
| bS® 
116 
117 
118@ 
119@ 
1200 
121@ 
122@ 
123@ 
124@ 
125@ 
126 
127 
128@ 


' XS=TIME INCREMENT BETWEEN POINTS DIV/S@ 
| YSS=VOLTS "BETWEEN POINTS DIV/25 

| VSS=VOLTS/DIV 

! TSSeaeeTRIGGER TIME SCALE 

\ 


Deltim=VAL(X$) 

Del vol=VAL(Ys$) 
Vscale=VAL( Vs ) 
Tscale=VAL(Ts$) 

| 

QUTPUT @Dfile;WS,Tscale,Vscale 
FOR IT=@ TO 1023 

Tam( I] )=I*Deltim 
Volt(I)=V(1)#Delvol 

NEXT I 

| 

PRINT “WRITING TO FILE" 
OUTPUT @Dfile;Tam(#) Volt(*#) 
ASSIGN @Dfile TO * 


BEEP 


| SHUHHHRERHHHHRAHHEHRAHUHAEGHHHHH HR et oat 


INPUT "DO YOU WANT ANOTHER WAVEFORM? Y=YES" , Ang 


IF An$="Y" 
| 

PRINT "PROGRAM TERMINATED" 
END 


THEN GOTO 430 
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1@ 

2 

3@ 

4@ 

5 

6@ 

7® 

BO 

3@ 

10 
119 
120 
13@ 
14@ 
15@ 
16@ 
17 
18@ 
19@ 
200 
21@ 
22® 
2 3@ 
240 
250 
Z6@ 
27 
280 
230 
30@ 
310 
32@ 
33@ 
340 
350 
36@ 
37 
380 
39@ 
400 
41@ 
472@ 
430 
440 
450 


APPENDIX B 
COMPUTER LISTING OF "PLOT" 


!'PROGRAM: PLOT 


THIS PROGRAM TAKES TWO FORMATTED BDAT FILES AND SUBTRACTS 

ONE FROM THE OTHER. IT WAS DESIGNED FOR USE USING WAVEFORMS 
OBTAINED ON THE TEXTRONIX 243@ DIGITAL OSCOPE. DATA IS TAKEN 
FROM THE FILES IN THE FOLLOWING ORDER. 


| 
! 
| 
| 
WS = FILE DESCRIPTION STRING 

| TSCALE = TIME SCALE OF THE SCOPE (SEC/DIV) 
VSCALE = VOLTAGE SCALE OF THE SCOPE (VOLIZiie 
! T(@-1023) = TIME ARRAY (SEC) 

V(@-10@23) = VOLTS ARRAY (VOLTS) 

| 

| 

| 

| 


THE SUBTRACTION RESULTS ARE THEN ABLED TO BE PLOTTED AS@Raee 
COPY ON THE THINK JET OR PEN PLOTTER. A STANDARO OPTIONS 
INCORPORATED FOR EASE WHEN USING THE PROGRAM FREQUENTLY WITH 
THE SAME PARAMETERS. THE PROGRAM ALLOWS FLEXIBILITY IN 
CHANGING THE PLOT’S PARAMTERS BEFORE THE HARD COPY. 


| SEES REE Ete 
INITIALIZATION 

| SETHE EEE SS Ett 
OPTION BASE @ 

DIM X_data( 1023) ,Y data(1@23),Y_axis_unitsS$[ 32] 

DIM Head_1b1${6@] ,Y_axis_1b1${60) ,X_axis_1b61$(32Z] 

DIM X_axis_units$[ 32] ,Clear_crt${Z] ,C${2] ,Sub head_l1b1 Sf 32) 
DIM X_char$[8] ,Vdata${15] ,WS[40]  VdataZ${[ 15) ,Y_chars[ 8] 

DIM X1(1023) ,Y1( 1023) ,X2( 1023) ,¥2( 1023) 
Clear_crt%$=CHRS( 255 )&CHRS( 75) 

C$=CHRS(255)8"K" 

| 


| SSH at st ee eet 
| 
GOSUB Load file 
INPUT "DO YOU WISH TO USE THE STANDARD PARAMETERS? Y=YES" ,Dec2$ 
IF Dec2$="Y" THEN 
GOSUB Stand params 
GOTO 42@ 
END IF 
GOSUB Initial_params 
AS="N" 
GOSUB Plot_data 
INPUT "WISH TO CHANGE CURRENT GRAPHICS PARAMETERS? Y=YES" ,DecZ$ 
IF Dec2$="Y" THEN GOSUB Change_params 
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AGO 
470 
480 
490 
500 
510 
520 
530 
540 
550 
560 
S70 
S80 
590 
500 
610 
620 
630 
640 
650 
560 
670 
680 
690 
700 
710 
720 
730 
T4Q 
75@ 
760 
77@ 
780 

790 
BOO 
BI? 
820 
830 
840 
850 
BE 
870 
880 
890 
900 
910 
920 
930 
940 
950 
360 


INPUT “ANOTHER DATA FILE PLOTTED ON SAME HARD COPY GRAPH?" ,Dectbs 
WHILE Dec6$="Y" 
GOSUB Load_file 
INPUT "LINE TYPE( 1=SOLID,4=DASH ,S=BIG DASH ,6=DASH-DOT” ,Line_num 
INPUT "PEN NUMBER 1 TO 8" ,Pen_num 
GOSUB Same_graph 
GOTO 460 
END WHILE 
INPUT "DO YOU WISH TO CHANGE ANY MORE PARAMTERS? Y=YES" ,Dec2$ 
IF Dec2$="Y" THEN GOTO 450 
INPUT "DO YOU WISH TO DUMP GRAPHICS? Y=YES" AS 
IF AS="Y" THEN GOSUB Plot_data 
INPUT "DO YOU WANT TO PLOT A DIFFERENT GRAPH? Y=YES" ,Dec7$ 
IF Dec7S="Y" THEN GOTO 35@ 
GOSUB End_program 
' 
| RHEAHRRARHIGINA HARRAH HHARHHHHHHH HHH Ht 
STANDARD PARAMETERS 


p38 S39:9-9°8:938201930502815 91058 1958381508: 65919191519+850:85518:5:8+9' 5103655585858 591919191868155818+55958-959551 810582583: 
! 


Stand params: |! 


Head lb1lS=W$S 

Sun _head_1b1$=Vdataés 
X_charS$e#"K" 

Y_char$="0D.D0D" 
Y_axis_I]lb1]$="Volts" 
X_axis_lb1¢="Relative Time" 
Y_axis_units$="Volts” 
X_axis_unitsS="microseconds 
Left=@ 

Right=!9#Tscale 
Bottom=-Vscalet4 
Top=Vscale#4 

X tic _sp=Iscale/72 
Y_tic_sp=Vscale/S 
X_ma)_ct=2 

Y_maj_ct=5 

Maj_tic=Z 

Xgoridtic= scale 
Ygridtic=Vscale 

Xor idma)}=1 

Yor idma)=1 

Majorid_sz2z#*Z 
Y_incmt<=2#Vscale 
X_incmt=4#Tscale 

Line_num= 1 

RETURN 

| 

Up 85 9391556359519 1818+ 9191510+858:0+018'9 2029081513 5595:950:0:9:0:8:0+01029:8:91655 191515181010: 85653:9:5:8:9:9:5:9:9:8:8:8: 8583; 
INITIAL PARAMETERS 

R19 399838:95910195919:9:5:915:019:101919191518:19'S:8:5:9: 9681918: 855'5:9:5: 915: 86:8:85 81581818: 955°51555°5:015:030:853:9:9:3. 5:5 


ot 
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97@ 


98@ Initial params: 


990 
1900 
1010 
102@ 
1030 
1042 
105@ 
1@6@ 
107@ 
1080 
1090 
119@ 
111@ 
112@ 
1130 
1142 
11S@ 
1160 
117@ 
1171 
118@ 
119@ 
12@0 
121@ 
122@ 
123@ 
1Z24@ 
125@ 
1260 
12°7@ 
1280 
1290 
1 300 
1301 
131@ 
132@ 
133@ 
1342 
135@ 
1360 
137 
138@ 
1390 
1400 
1410 
142@ 


X tic_so=@ 

Y-t1¢-so—¢ 

X maj _ct=!] 

Y_maj_ct=] 

May tic=2 

Xor idtic=@ 

Yoridtic=@ 

Xaor 1dmaj=1 

Yor idma)=! 

Majygrid_s2=2 

INPUT “GRAPH HEADING LABEL ,32 CHAR. MAX” Head_!b1$ 

INPUT “GRAPH SUB LABEL ,32 CHAR. MAX" ,Sub_head_1bl1$% 

INPUT "IMAGE FOR X AXIS NUMERICS, 8 CHAR. MAX" |_X_chars 
INPUT “IMAGE FCF Y AXIS NUMERICS, 8 CHAR. MAX" ,Y_chars$ 
INPUT "Y AXIS LrBEL,32 CHAR. MAX”, Y_axis_l1bl1$ 

INPUT "X AXIS LASEL,32 CHAR. MAX" ,X_axis_1bI]$ 

INPUT "“Y AXIS UNITS,32 CHAR. MAX" Y_axis_units$ 

INFUT “X AXIS UNITS,32 CHAR. MAX" |X _axis_units$ 

DIsP “WINDOW LOOKED AT,LEFT,RIGHT ,BOTTOM, TOP”; 

INPUT Left,Right,Bottom, Top 

INPUT "X AXIS TICK SPACING, DEF=O(NO TICKS)” |X _tic_sp 
INPUT “Y AXIS TICK SPACING ,DEF=O(NO TICKS)” ,Y_tic_sp 
INPUT “NUM TIC INTVLS BIWN MAJ TICS X-DIR,DEF=1" ,X_maj_ct 
INPUT “NUM TIC INTVLS BTWN MAJ TICS Y-DIR,DEF=1",Y_maj_ct 
INPUT “SIZE OF MAJOR TICKS, DEF=2" ,Maj_tic 

INPUT "GRID X-TIC SPACING,DEF=@(NO TICKS)”, Xgridtic 

INPUT “GRID Y-TIC SPACING, DEF=@(NO TICKS)" ,Yoridtic 

INPUT “NUM TICKS BIWN MAJOR TICKS Y-DIR, DEF=1" , Ygridmaj 
INPUT "NUM TICKS BIWN MAJOR TICKS X-DIR, DEF=1" , Xgridmaj 
INPUT "SIZE OF MAJOR GRID TICKS,DEF=2" ,Majgorid_sz 

INPUT "Y-AXIS INCREMENT” ,Y_incmt 

INPUT "X-AXIS INCREMENT” |X _incmt 

DISP “LINE TYPE( 1=SOLID,4=DASH,5=BIG DASH,6=DASH-DOT)"; 
INPUT Line_num 

| 

RE TURN 

| 

B18: 0:950:850:510:8:0:0:0:5+8:0:8:8°9:9:2:8:9) 8:85 
! LOAD DATA FILE 


SBoeStSr: Sr PS et He St Se Sr Sr: OreGr: Sr  HrcSi Gre Or GS: Ge,Hr1 S22 OH2Or2Gr Gi A osorgi Ori S1G:Sreraieace 








Load file: |! 


: | 
INPUT “INPUT NAME OF DATA FILE 1” ,Vdata$ 
INPUT “INPUT NAME OF DATA FILE 2 “ ,VdataZz$ 
ASSIGN @Df1 TO Vdata$ 


WZ 


1430 
1449 
145@ 
1469 
147 
148@ 
149@ 
1500 
151@ 
1520 
1530 
154@ 
1550 
1560 
1572 
1580 
1590 
160@ 
1610 
162@ 
163 
164@ 
1650 
1669 
167 
168Q 
16390 
1700 
1710 
17270 
1732 
174Q 
175@ 
176 
1770 
1780 
1739@ 
180d 
181@ 
182@ 
183@ 
184@ 
1859 
1860 
187@ 
1880 
1890 
1300 


ENTER @DfF13WS$,Tscale,Vscale,XI(*#) YIC#) 
ASSIGN @Dfl TO * 

ASSIGN ®DfZ2 TO Vdata7$ 

ENTER @Df23WS,Tscale,Vscale,X2(*) ,Y2(*#) 
ASSIGN @Df2 TO #* 

MAT X_data= Xl] 


MAT Y_data= Y1-Y2 
PRINT "DATA FILES READ" 
BEEP 

RETURN 


I 

UB 825:9:9:550:9:959:8:9:81559:3:0:0:9)008:8: 8-525 :8:518:5:950:5:601819191018:8:8181875:015:8:985959:515:9:9:8:858:5:5:5:5.8 385: 
| PLOT DATA 

UB + 23029305999186018369819995910+65918585018192950203995°5:8151598:5291959:919:9751616:555:5:5+9+536-9:59019.5:9: S1s3S.s2508 258s. 
1 


Plot_data: | 


IF AS="Y" THEN 
OUTPUL KBO;Clear_crts; 
GINIT 
medvterels CRI. “INTERNAL” 
GOTO 1772 
eno LE 
itmel is THIS 16 BEA HARDCOPY PLOT? Y=YES”" ,DecI$ 
IF Dec!$="Y" THEN 
GINIT 
PLOTTER IS 705," HPGL” 
GOTO 1772 
ESE 
OUTPUT KBD;:Clear_crt$; 
GINIT 
PLOTTER IS CRT," INTERNAL" 
GOTO 1772 
END IF 
GRAPHICS ON 
PEN 1 
X_Qdu_max=100#MAX( 1 RATIO) 
Y_gdu_max=100#MAX(1,1/RATIO) 
LORG 5 
CSIZE 4 
MOVE .48*xX gdu_max, .93#Y_Qdu_max 
LABEL Head_lbl$ 
MOVE .48*xX_gdu_max,.89#Y_qdu_max 
esac. 3 
LABEL Sub head_lbl$ 
DEG 
LDOIR 9@ 
CSIZE 4 
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191@ 
192@ 
1930 
194@ 
1950 
196@ 
1972 
198@ 
1990 
2 OOO 
2010 
2OZ@ 
2030 
2040 
2050 
2260 
2O7@ 
298O 
2090 
2100 
2110 
Z212@ 
2130 
214@ 
2150 
216¢ 
217@ 
2180 
219@ 
220@ 
2219 
2220 
2232 
224@ 
2252 
2260 
2270 
2280 
2290 
Z 300 
231® 
2320 
2330 
2340 


MOVE .@4#xX gdu_max,.5S5*#Y_gdu_max 
LABEL Y_axis_I1lb1$ 
LORG 5S 
LDIR @ 
MOVE .48#X gdu_max,.15#Y_ gqdu_max 
CSIZE 4 
LABEL X_axis_1lbl1$% 
MOVE .48#X gdu_max,.19#Y gdu_max 
CSIZE 3 
LABEL “("3;X_axis_units#;")" 
MOVE .70#X_gdu_max,.74#Y_ gdu_max 
VIEWPORT .15*X_gdu_max,.85#X_gdu_max ,.25#Y_gdu_max ,.86#Y_gdu_max 
WINDOW Left ,Right ,Bottom, Top 
AXES X_tic_sp,Y_tic_sp,Left,Bottom,X_maj}_ct,Y_maj_ct,Ma)_tic 
AXES X_tic_sp,Y_tic_sp,Right , Top,X_maj_ct,Y_maj_ct ,Maj_tic 
IGRID Xgridtic, Yoridtic,Left ,Bottom, Xgridmaj ,Yoqridmaj ,Majgrid_sz 
COVE OFF 
CSIPZE 3,.75 
LORG 6 
FOR T=Left TO Right STEP X_incmt!LABELING OF X AXIS 
MOVE I ,(Bottom-.@4+#( Top)) 
LABEL USING X_charS; l74ere-6 
NEXT I 
LORG 8 
FOR I=Bottom TO Top STEP Y_incmt!LABELING OF Y AXIS 
MOVE (Left-.4#Tscale) I 
LABEL USING Y_char$;I 
NEXT I 
PENUP 
PEN 2 
LINE TYPE Line_num 
CLIP ON 
FOR I=@ 10 1023!PLOTTING OF DATA 
PLOT X_data(I),Y_data(!) 
NEXT I 
PE NUP 
CLiP OFF 
IF AS="Y" THEN OUMP GRAPHICS #701 
RE TURN 
j 
UB 10:5: 9:9:8-9:9:5:9:9:9:0'5+950:9:810:89:8:8: 9+) 
l SAME GRAPH 







2350 Same_graph: |! 


2 360 
2370 
2380 


PEN Pen_nunm 
LINE TYPE Line_num 
CLIP ON 
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2390 
2400 
2410 
2420 
2430 
244 
2450 
2460 
2470 
2480 
2490 
2500 
251 
2520 
2530 
2540 
2550 
2560 
2570 
2580 
2590 
2600 
261 
2620 
2630 
2640 
2652 
2660 
Z267@ 
2680 
(2690 
2702 
271@ 
2722 
2730 
2742 
2750 
2760 
2770 
2780 
27292 
2 BOW 
2812 
2820 
2830 
2842 
2850 
2 860 


FOR T=O@ [FO 1023 
PLOT X_data(I),Y data(!I) 
NEXT I 
CLEP DFE 
PENUP 
RETURN 
i 
UB 3'9:9:6:9191838:9:098:9:6538:5: 91558: 93859:978:518151810:9:9:9:5:9'S:558)955°8:8:91818:9:9:5:8:3:9:5:9:8:8:515:5:9:91928:9: 3): 
CHANGE GRAPHICS PARAMETERS 
| HHHHERHRHHHRHEIIR HR a Ee at 
| 
Change_params: ! 
OUTPUT KBO;Clear_crt$; 
PRINT “SELECT FROM THE FOLLOWING TO CHANGE” 
PRINT "CASE 1: GRAPH HEADING" 
PRINT "CASE 2: Y AXIS LABEL" 
Piaend CASE 325X AXIS LABEL” 
PRINT "CASE 4: Y AXIS UNITS” 
PRINT “CASE S&S: X AXIS UNITS" 
PRINT "CASE 6: WINDOW UNITS" 
bien. CAGE 7: X AXIS TICK SPACING” 
een CAGE 8: Y AXIS TICK SPACING" 
mein) CASE Se TICK LNTVin BIWN MAJOR TICKS IN Xe OIRECTION” 
PRINT “CASE i@: TICK INTVL BTWN MAJOR TICKS IN Y OIRECTION" 
Pree Ghse Tl: SkZE OF “MAJOR” TICKS” 
PRINT meGASE 17: XK @ICK GREDeSPACING’ 
PRINT “CASE 13: °¥° TICK GRID SPACING” 
PRINT "CASE 14: GRID TICK INTVL BTWN MAJOR TICKS IN X OIRECTION" 
P&INT "CASE 15: GRID TICK INTVL BTWN MAJOR TICKS IN Y DIRECTION" 
PRINT “CASE 16: MAJOR GRID TICK SIZE" 
PRINT “CASE 17: Y INCREMENT" 
PRINT "CASE 18: X INCREMENT" 
PRINT “CASE 19: LINE TYPE” 
PRINT “CASE 20: GRAPH HEADING SUB LABEL" 
PRINT “CASE 21: X AXIS NUMERIC IMAGE" 
PRINT "CASE 22: Y AXIS NUMERIC IMAGE" 
| 
INPUT "DO YOU DESIRE TO CHANGE ANY PARAMETERS? Y=YES" ,Dec4S$ 
IF Dec4s="Y" THEN 
INPUT "ENTER PARAMETER NUMBER" ,Change_num 
ELSE 
GOTO 3560 
END IF 
SELECT Change_num 
CASE | 
PRINT "CURRENT HEADING LABEL" ,Head_1b1$ 
INPUT "NEW HEADING LABEL" ,Head_1b1$ 
GASE 2 


S 


287@ 
2880 
2830 
290 
291@ 
292@ 
293 3@ 
294@ 
2950 
2960 
29°7@ 
298@ 
2390 
3000 
3@1@ 
3022 
3030 
3@40 
3050 
3Q6® 
3@7@ 
3080 
3090 
3100 
311@ 
oe 
3130 
3148 
Stade 
3160 
317@ 
318@ 
313@ 
3200 
3210 
322@ 
32 3@ 
324@ 
3250 
326@ 
327@ 
3280 
329@ 
3300 
331@ 
332@ 
333@ 
334@ 


PRINT 
INPUT 
Cassese 
PRINT 
INPUT 
CASE 4 
PRINT 
INPUT 
CASE S 
PRINT 
INPUT 
CASE 6 
PRINT 
INPUT 
COSEma 
PRINT 
INPUT 
CASE 8 
PRINT 
INPUT 
CASE <3 
PRINT 
INPUT 
CASE 10 
PRINT 
INPUT 
Cnsta 
PRINT 
INPUT 
CASE 12 
PRINT 
INPUT 
CASE Ps 
PRINT 
INPUT 
CASE 14 
PRINT 
INPUT 
CASE 1S 
PRINT 
INPUT 
CASE 16 
PRINT 
INPUT 
CASE 17 
PRINT 
INPUT 
CASE 18 


“CURRENT Y AXIS LABEL ,Y¥caxneu ils 
“NEW Y AXIS LABEL” ,Y_axis_1b1$% 


“CURRENT X AXIS LABEL" ,X_axis_1bl$ 
“NEW X AXIS LABEL" ,X axis_1bl1$ 


“CURRENT Y AXIS UNITS" Y_axis_units$ 
"NEW Y AXIS UNITS"  Y_axis_units$ 


“CURRENT X AXIS UNITS” X_axis_units$ 
"NEW X AXIS UNITS" X_axis_units$ 


“CURRENT WINDOW" Left Right ,Bottom,Top 
"NEW WINDOW" Left, Right ,Bottom, Top 


“CURRENT X AXIS TIC SPACING" ,X_tic_sp 
“NEW X AXIS TIC SPACING" ,X_tic_sp 


“CURRENT Y AXIS TICK "SPACING | Y ticesp 
“NEW Y AXIS TIC SPRGING =a vie sso 


“CURRENT TIC INTVL BIWN MAJ TICS X DIR" ,X_maj_ct 
“NEW TICK INTVL BIWN MAJOR TICKS X DIR" ,X_maj_ct 


“CURRENT TICK INTVL BTWN MAJOR TICKS Y DIR", Y_maj_ct 
“NEW TICK INTVL BTWN MAJOR TICKS Y DIR", Y_maj_ct 


“CURRENT MAJOR TICK SIZE" ,Maj_tic 
"NEW MAJOR TICK SIZE" ,Maj_tic 


“CURRENT X TICK GRIO SPACING" ,X_grid_tic 
“NEW X TICK GRIO SPACING", X_grid_tic 


“CURRENT Y TICK GRID SPACING” ,Y_grid_tic 
"NEW Y TICK GRID SPACING" ,Y_grid_tic 


"CURRENT GRID TIC INTVL BTWN MAJ TICS X DIR" ,X_grid_ma}j 
"NEW GRID TIC INTVL BTWN MAJ TIC X DIR" ,X_grid_maj 


“CURRENT GRID TIC INTVL BTWN MAJ TICS Y DIR" ,Y_grid_maj 
“NEW GRID TIC INTVL BTWN MAJ TIC Y OIR",Y_grid_maj 


“CURRENT MAJOR TICK SIZE” ,Maj_grid_sz 
“NEW MAJOR TICK SIZE" ,Maj_orid_sz 


“CURRENT Y INCREMENT”, Y_incmt 
“NEW Y INCREMENT” ,Y_incemt 
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3350 
3360 
BRA’) 
3380 
3390 
S391 

3400 
341@ 
3420 
3430 
3440 
3450 
3460 
347@ 
3480 
3490 
3520 
3519 
Sal 
eB, 
354@ 
3550 
355@ 
3570 
358 
353@ 
3600 
3610 
3620 
3630 
3642 
365@ 


PRINT “CURRENT X INCREMENT" , X_incmt 
PRINT "NEW X INCREMENT", X_incmt 
(Aor sed 
PRINT “CURRENT LINE TYPE" ,Line_num 
DISP “NEW LINE TYPE(1=SOLID,4=DASH ,S=BIG DASH,6=DASH- DOT"; 
INPUT Line_num 
CASE 20 
PRINT “CURRENT HEADING SUB LABEL” ,Sub_head_1b1$ 
INPUT “NEW HEADING SUB LABEL" ,Sub_head_1lbl$ 
GAse 721 
PRINT “CURRENT X AXIS NUMERIC IMAGE” ,X_char$ 
INPUT “NEW X AXIS NUMERIC IMAGE" ,X_char$ 
CASE 22 
PRINT “CURRENT Y AXIS NUMERIC IMAGE” ,Y_char$ 
INPUT "NEW Y AXIS NUMERIC IMAGE” ,Y_char$ 
END SELECT 
INPUT "DO YOU WISH TO CHANGE ANOTHER PARAMETER? Y=YES" ,DecS$ 
IF DecS$="Y" THEN 
GOTO 2510 
bluGE 
GOSUB Plot_data 
END IF 
RETURN 
{ 
I HHRUHRERHHHUHHRATHEARAHRNHHHERAE HAR ERER HHA ARH SASH 
1 ENO THE PROGRAM 
| SHHRRURHHAH GR RHRHRERARA NARHA HRA SSH 
' 
End_program: |! 
GRAPHICS OFF 
PRINT “END OF PROGRAM" 
END 
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